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ABSTRACT 

We investigate the effects of varying the intensity of the primordial magnetic seed 
field on the global properties of the galaxy population in ideal magnetohydrodynamic 
cosmological simulations performed with the moving-mesh code AREPO. We vary the 
seed field in our calculations in a range of values still compatible with the current 
cosmological upper limits. We show that above a critical intensity of ~ 10“® G, the 
additional pressure arising from the field strongly affects the evolution of gaseous 
structures, leading to a suppression of the cosmic star formation history, which is 
stronger for larger seed fields. This directly reflects into a lower total galaxy count 
above a fixed stellar mass threshold at all redshifts, and a lower galaxy number density 
at fixed stellar mass and a less massive stellar component at fixed virial mass at all 
mass scales. These signatures may be used, in addition to the existing methods, to 
derive tighter constraints on primordial magnetic seed field intensities. 
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1 INTRODUCTION 

Much progress has been made to understand the evolution 
of magnetic fields in the cosmological context, and espe¬ 
cially on investigating mechanisms through which B fields 
are set to their present-day strength in galaxies and galaxy 
clusters. Although there are strong indications that dynamo 
processes play a crucial role in amplifying pre-existing B 
fields (Hanasz et al. 2009; Schober et al. 2013), the origin 
of these seed fields still remains uncertain. In particular, it 
is unclear whether these seeds are of cosmological origin, 
i.e. produced at shocks (Ryu et al. 1998) or during infla¬ 
tion and phase transitions in the early Universe (Widrow 
et al. 2012), or generated by stars in proto-galaxies (Schle¬ 
icher et al. 2010) and subsequently ejected by galactic winds 
(Volk & Atoyan 2000) and active galactic nucleus activity 
(Furlanetto & Loeb 2001). 

It is intrinsically difficult to observe magnetic fields out¬ 
side haloes or at high redshift with the current generation of 
instruments (Beck 2007). The detection of magnetic fields 
outside collapsed structures is particularly challenging be¬ 
cause the expected intensities ( ^ 1 nG) are well below the 
ones observed in galaxies and galaxy clusters (~ lfj,G). 
Hence, most of the constraints on primordial magnetic fields 
are upper limits. These are in turn derived by considering 
how the predicted effect of such seed fields would compare to 
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actual observations of more well-determined quantities such 
as the cosmic microwave background power spectrum or po¬ 
larization (Barrow et al. 1997; Durrer et al. 2000; Jedamzik 
et al. 2000; Kahniashvili et al. 2009, 2010; Planck Collabora¬ 
tion XIX 2015), the isotropy in the distribution of ultrahigh 
energy cosmic rays (e.g. Lee et al. 1995), Big bang nucleosyn¬ 
thesis abundances (Grasso & Rubinstein 1995, 2001), Farda- 
day rotation measure, Ly a forest and Sunyaev-Zel’dovich 
effect statistics (Blasi et al. 1999; Shaw & Lewis 2012). In¬ 
terestingly, most of the limits converge towards a maximum 
present-day field strength of ~ 1 — 10 nG. Moreover, possi¬ 
ble lower limits ( > 10“^® G) have been inferred from y-ray 
and TeV blazars observations (Neronov & Semikoz 2009; 
Neronov & Vovk 2010). 

In this Letter, we examine how strong the effects of cos¬ 
mological B fields (i.e. fields generated in the early Universe) 
on the global properties of the galaxy population are as a 
function of the field strength, and determine what is the crit¬ 
ical intensity at which they become noticeable. We run a set 
of cosmological magnetohydrodynamic (MHD) simulations 
with varying seed field intensities and study the changes 
in the cosmic star formation density, stellar mass function, 
stellar-to-halo mass relation, the evolution of the mean mag¬ 
netic field and total galaxy number. The Letter is organized 
as follows. In Section 2, we present analytic considerations 
to relate the seed field used in our runs to the current up¬ 
per limits and we estimate at what intensity those fields are 
expected to significantly impact the dynamics of the gas. 
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We introduce the setup of our simulations in Section 3 and 
we describe the main results of our analysis in Section 4. 
Finally, we draw our conclusions in Section 5. 


2 ANALYTIC CONSIDERATIONS 


In our simulations, the initial seed magnetic field is uniform 
and threads the computational domain along a predefined 
direction (see also Sect. 3). However, in the majority of ear¬ 
lier work, upper limits for cosmological magnetic fields are 
inferred by assuming that those fields are described by a 
power spectrum of the form (see e.g. Kahniashvili et al. 
2009) 

PB{k) = ABk"^ k<kn, (1) 


where Ab is a normalization constant, ub is the spectral in¬ 
dex and ku is the so-called Alfven damping scale (Jedamzik 
et al. 1998) above which it is assumed that no power is 
present. Therefore, it is important to establish a direct con¬ 
nection between the uniform field intensities used in this 
work and those derived through equation (1). Before pro¬ 
ceeding, we note that the value of ub is uncertain (usually 
—3 < riB SS 2; Kahniashvili et al. 2010) and can have a 
significant impact on the value of the inferred intensities 
of cosmological B fields. It is also to avoid this additional 
source of uncertainty that we opted for the simpler uniform 
field setup in our runs. 

The information about the intensity of cosmological 
seed B fields is contained in the normalization constant Ab, 
which is linked to a smoothed comoving field amplitude on 
a comoving scale A via a convolution with a Gaussian kernel 
through (Planck Collaboration XIX 2015) 
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where P is Euler’s gamma function. From equation (2) it is 
easy to connect the smoothed field amplitude at two arbi¬ 
trary scales Al and A 2 as 


Bx 2 = Bxi 



(3) 


In our simulation setup, the initial seed field is coherent 
on a 25h~^Mpc ~ 37 Mpc scale. Hence, a smoothed field 
Bi Mpc = 1 nG - a representative value of current upper 
limits (Planck Gollaboration XIX 2015) - translates into 
BsTMpc — 1.2 X 10“"^ nG for ub = 2, and Bstmpc — 0.84nG 
for riB = —2.9. 

A more physically motivated way to express the inten¬ 
sity of a magnetic seed field whose power spectrum satisfies 
equation (1) is to replace Bx with the effective field contain¬ 
ing the same magnetic energy density as the one obtained by 
integrating the power spectrum over all scales. In this way, 
the information at all spatial scales is used and the resulting 
field intensity can be compared in a more meaningful way 
to the uniform seed field that we have adopted in our calcu¬ 
lations. The effective field is defined as (Kahniashvili et al. 
2010, equation 5) 


RA ( fcDA )("'^+^)/2 

ri/2(„B/2-H5/2)' 


(4) 


By applying equation (4) in the case of the ub interval that 
we have examined above, we obtain an effective field Beff — 
1.22 nG and Bes ~ 240 nG for ub = —2.9 and ub = 2, 
respectively. These values have been derived for BiMpc = 
1 nG and by adopting the expression for ko given in equation 
(3) of Planck Gollaboration XIX (2015). The maximum seed 
field value that we employ in our simulations is consistent 
with those limits. 

As we mentioned in the Introduction, the majority of 
the currently available upper limits translate into effective 
seed fields with intensities of the order of 1 — 10 nG. We 
now present two simple arguments showing that significant 
effects on the evolution of the galaxy population must be 
expected if the Universe is permeated by those (or stronger) 
fields. We start by noting that in ideal MHD, the magnetic 
field flux is conserved. This leads to the scaling relation 
(Grasso & Rubinstein 2001) 

where Bo is the intensity of the initial seed field and 5 is 
the density contrast between the average gas density within 
galaxies and that in the intergalactic space. Equation (5) 
states that the final field intensity observed nowadays in 
galaxies (~ a few /rG) can be entirely accounted for by adi¬ 
abatic compression of an initial seed field of 1 nG. However, 
several numerical studies (see Marinacci et al. 2015, and 
references therein) have shown that additional processes, 
such as small-scale dynamo and turbulent or shearing flows, 
play a key role in amplifying the magnetic field to its final 
strength. It is thus plausible to assume that seed fields of 
~ 1 nG would be over-amplified by the presence of these 
(small-scale) mechanisms. This over-amplification would in 
turn leave its imprint on galactic properties by enhancing 
the dynamical role of magnetic fields. 

Magnetic fields start to play a role in the gas dynamics 
if the magnetic pressure becomes comparable to the thermal 
gas pressure, viz. 

||B||">8^p-, (6) 

fJjTTlp 

where ||B|| is the magnetic field intensity p and T the gas 
density and temperature, k the Boltzmann constant, p the 
gas mean molecular weight and rrip the proton mass. By 
considering a fully-ionized gas of primordial composition 
{X = 0.76), we can rewrite equation (6) in terms of the 
comoving magnetic field intensity as 

where h = Bo/lOOkm s“^ Mpc“^, fib is the baryon density 
and « the redshift. 

Equation (7) shows that the transition to magnetically 
dominated gas flows occurs at a critical seed magnetic field 
strength between 1 and 10 nG, in line with the estimates 
and the observational constraints discussed earlier. In par¬ 
ticular, in our runs the gas reaches temperatures of ~ 10'* K 
at the end of the reionization, which occurs at Zreion = 6 for 
our choice of the UV background (Faucher-Giguere et al. 
2009). This implies a critical seed B field IjiJjj ^ 1.44 nG. 
Therefore, we anticipate a significant effect on galactic evo¬ 
lutionary processes at these B field intensities. 
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Figure 1. Left: evolution of the volume-weighted rms B field strength. Above the critical B field value of 10“® G the field intensity is 
set by flux conservation and little turbulent amplification is present. Centre: evolution of the cosmic star formation rate density (SFRD) 
together with a compilation of observational data taken from Behroozi et al. (2013). For larger seed magnetic field strengths there is an 
increasing suppression of the SFRD, which is present at all times for seed fields > 10“® G. Right: evolution of the total galaxy number 
above a stellar mass threshold of 10®Mq. The decrease in the total number of objects for seed field > 10“® G is visible at all times. 


3 SIMULATIONS SETUP 


Our simulations evolve a uniformly sampled box from the 
initial redshift z = 127 to the present day. We adopt a A 
cold dark matter cosmology with parameters flm = 0.302, 
fib = 0.04751, flA = 0.698, erg = 0.817, = 0.9671, 

and Ho = 68 km s“^Mpc“^ derived by an independent re¬ 
analysis of Planck data by Spergel et al. (2015). All the runs 
are performed with the moving-mesh code AREPO (Springel 
2010), together with a module for ideal MHD (Pakmor & 
Springel 2013) based on the 8-wave Powell et al. (1999) for¬ 
malism to control divergence errors. The code is comple¬ 
mented with a galaxy formation physics model developed 
for the ILLUSTRIS simulation suite (Vogelsberger et al. 2014), 
which includes as major components: (1) primordial and 
metal-line cooling with self-shielding corrections (Rahmati 
et al. 2013), (2) an updated version of the Springel & Hern- 
quist (2003) model for star formation with a Chabrier (2003) 
initial mass function, (3) stellar evolution, gas recycling and 
chemical enrichment, (4) a kinetic galactic wind scheme in 
which the wind speed is determined by the local dark matter 
velocity dispersion in haloes (see e.g. Puchwein & Springel 
2013), (5) black hole seeding, accretion, merging and feed¬ 
back (see also Springel et al. 2005). In this work we use the 
fiducial settings of the model presented in Vogelsberger et al. 
(2013). 

We sample a 25 h“^Mpc box size with 2 x 256^ reso¬ 
lution elements. Magnetic fields are added at the beginning 
of the calculation as a uniform seed field throughout the 
box. We choose this initial setup because it is the simplest 
divergence-free configuration to implement and has also the 
advantage that the field is specified once its intensity is fixed 
(but see Marinacci et al. 2015, for a discussion on how the 
initial field direction affects its properties). We keep the di¬ 
rection of the field fixed in all simulations along the z-axis, 
but change the (comoving) strength of the initial seed in the 
interval [10~^^ — 10“®] G. 


4 RESULTS 

The left-hand panel of Fig. 1 presents the time evolution of 
the (volume-weighted) rms B field intensity as a function of 
redshift. The initial seed field intensity is indicated in the 
legend. At low seed fields (< 10“®G), after a first phase 
dominated by the expansion of the Universe, in which the B 
field strength declines oc (1 -I- zfi (grey dashed line), there 
is an upturn of the mean B field intensity at z ~ 6, after 
which the field is amplified by structure formation. The am¬ 
plification saturates quickly, leaving behind an average field 
of ~ 10“^ fiG, whose amplitude slowly declines with time. 
The situation is different for seed fields ^ 10“® G. In those 
cases, the amplification of the field due to build-up of cos¬ 
mic structures is less prominent. Indeed, the magnetic field 
strength steadily drops with time and adiabatic compres¬ 
sion alone is largely sufficient to set the B field intensities 
observed in the present-day Universe (see eq. [5]). 

The increased magnetic field strength for large seed 
fields implies a larger magnetic pressure acting on the gas, 
which in turn provides more support against gravitational 
collapse, potentially affecting the efficiency of its conversion 
into stars. We illustrate this in the central panel of Fig. 1 
in which we show the cosmic star formation rate density 
(SFRD) as a function of redshift compared to a compila¬ 
tion of observational determinations taken from Behroozi 
et al. (2013). We can roughly divide the simulation set in 
two groups, depending on whether the initial seed field is 
below or above the critical value of 10“® G, which is very 
close to the one that we have determined analytically. Be¬ 
low this value, the magnetic field is still dynamically negli¬ 
gible (see eq. [7]) and the effects on the SFRD are minimal. 
Conversely, seed fields 10“® G have a significant impact 
on gas dynamics. The result is a noticeable reduction of the 
global SFRD, which is more prominent the larger the seed 
field and increasingly falls short of the observational findings 
especially at the star formation peak at 2 ~ 2, while recon¬ 
ciling the simulations to the observed level of star formation 
at 2 ~ 0. The simulation with seed field 10“® G exhibits an 
intermediate behaviour. At 2 ^ 3, it follows very closely the 
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same trend of the simulations with a lower seed field. A re¬ 
duction in the SFRD (about a factor of 2 with respect to 
the runs with seed field < 10“® G) only occurs after 2 ; ~ 3. 

The right-hand panel of Fig. 1 displays the evolution 
of the total number of galaxy above a stellar mass thresh¬ 
old of 10®Mq. The stellar mass is determined as the mass 
in stars contained within twice the stellar half-mass radius 
of each object. From the plot, it can be readily appreciated 
that the number of galaxies grows very rapidly for 2^6 and 
by 2 ~ 2 (approximately at the peak of the cosmic SFRD) 
already a significant fraction ( 60%) of the total objects 

has formed. The suppression of the cosmic SFRD due to in¬ 
creasing B seed fields intensities translates into a decrease 
of the total number of galaxy formed. The decrease is more 
pronounced for larger B field strengths and particularly no¬ 
ticeable for all the runs with fields above the critical inten¬ 
sities value of 10~® G. In those cases, the deficiency in the 
total galaxy number relative to simulations with seed field 
strengths below the critical value - which closely track one 
another - is present at all redshifts. Again, the 10“® G sim¬ 
ulation shows an intermediate behaviour starting to slightly 
deviate from the trend followed by lower seed field runs only 
after 2 ~ 2. We explicitly checked that the general trends 
discussed above are roughly independent from the partic¬ 
ular value of the stellar mass cut adopted, at least up to 
10^° Mq where minor differences in the final total number 
of objects formed start to be present also in runs with seed 
field intensities < 10“® G. 

The top panel of Fig. 2 shows the redshift zero galaxy 
stellar mass function for all simulations together with ob¬ 
servational values inferred from SDSS data (Bernardi et al. 
2013). Due to the inability of small haloes to accrete highly 
magnetized gas, we expect a global reduction of the num¬ 
ber of galaxy per unit stellar mass for increasing seed field 
intensities. The plot fully confirms these expectations and 
shows that the effect is stronger the larger the initial seed 
field. For seed fields ^ 10“® G, the mass functions agree, 
within the scatter, over all the considered mass range. For 
seed fields above 10“® G, the reduction in number density 
is clearly visible up to M* ~ Mq, where the number 

of galaxies is underpredicted compared to the observations. 
At a stellar mass of ~ 10^° ® Mq, the discrepancy among 
the different runs is less prominent, an effect likely due to 
the fact that these haloes are massive enough to power their 
star formation even in the presence of a high B field. For 
larger M* (up to ~ 10^^ ® Mq) the differences grow again, 
signalling that the building blocks of these large structures, 
which accrete a significant fraction of their mass through 
mergers, contain on average less stars for large seed fields. 

In the bottom panel of Fig. 2, we present the redshift 
zero stellar-to-halo mass relation as a function of the virial 
mass. For comparison, we also show the same relation de¬ 
rived via abundance matching (Behroozi et al. 2013). Again, 
the net effect of an increase of the seed field strength is a 
reduction of the ratio between stellar and halo mass over 
the whole range of virial masses in the simulations. The 
general division in two groups discussed previously can also 
be applied here. In particular, simulations with seed fields 
< 10“® G track each other quite well, while for seed fields 
^ 10“® G, the stellar-to-halo mass ratio becomes smaller for 
increasing value of the initial seed magnetic field. At low halo 
masses, the conversion efficiency of gas into stars is consider- 
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Figure 2. Top: galaxy stellar mass function compared to SDSS 
observations (Bernardi et al. 2013). The suppression of star for¬ 
mation causes a decrease of the number of objects at all mass 
scales. The effect is stronger for larger seed field intensities. Bot¬ 
tom: stellar-to-halo mass relation contrasted to the abundance 
matching predictions of Behroozi et al. (2013). At any given virial 
mass, the stellar content of haloes decreases for increasing seed B 
field strengths. 

ably reduced for the highest seed magnetic fields, for which 
virial masses ^ 10^ Mq are required on average for a halo 
to host stars. The maximum conversion efficiency is found 
for virial masses ~ 10 ^^ ® Mq, about ~ 0.5 dex larger than 
the analogous relation found in illustris and inferred from 
abundance matching. Moreover, the peak of the relation ap¬ 
pears to be shifting, albeit only slightly, towards large halo 
masses. At the high-mass end (M2oo,crit ~ 10 ^®'® Mq), the 
discrepancy among curves gets smaller, but the main trend 
with the seed magnetic field intensity is still present. 


5 SUMMARY AND CONCLUSIONS 

We have investigated how the presence of cosmological mag¬ 
netic fields affects the global properties of the galaxy popula- 
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tion using cosmological MHD simulations. Our main results 
can be summarized as follows. 

(i) For seed B fields f 10“® G, the increased magnetic 
pressure hampers gas accretion into low-mass dark matter 
haloes. This causes a reduction of the cosmic SFRD, which 
becomes more marked for increasing seed field intensities. 

(ii) The suppression of the cosmic SFRD also results in a 
reduction of the total galaxy number at all redshifts and of 
present-day galaxy stellar mass function in all stellar mass 
bins and in particular for M^, f 10^° Mq . Again, this effect 
is stronger for larger seed B fields. 

(iii) At fixed virial mass, haloes tend to contain signifi¬ 
cantly less stars for increasing seed field intensities. This is 
particularly evident for the less massive haloes. Indeed, for 
the largest seed field considered in the simulations, virial 
masses f Mq are required for a halo to host stars. 

Our simulations neglect any change that the presence of 
a seed magnetic field can induce on the initial matter power 
spectrum. Several authors (e.g. Sethi & Subramanian 2005, 
2009; Yamazaki et al. 2006; Fedeli & Moscardini 2012) have 
argued that seed magnetic fields can affect the distribution 
of power especially at small/intermediate scales, at which 
the clustering of the matter would be enhanced. By not con¬ 
sidering these variations to the matter power spectrum, we 
might overestimate the effect of large B seed field on small 
haloes in our runs. Moreover, our cosmological volume is 
fully magnetized at the initial redshift, which again may re¬ 
inforce the effects of B fields on the gas with respect to a 
scenario in which the field is first generated within galaxies 
and then ejected in the intergalactic space. 

Another point to bear in mind is that degree of reduc¬ 
tion of the stellar content and star formation in haloes might 
also depend on the choice of the feedback model adopted 
in the simulations. The fiducial illustris feedback model 
used in our runs is at the upper end of what is energetically 
plausible, in particular for the generation of galactic winds. 
Therefore, it is conceivable that (strong) seed magnetic fields 
would cause a smaller suppression of the star formation rate 
and galaxy stellar masses with less efficient feedback loops, 
helping to alleviate some of the tensions with the observa¬ 
tions that we have discussed. 

Even with these caveats, from our runs emerges a con¬ 
sistent picture in which primordial magnetic seed fields and 
galaxy properties are closely related. In particular, primor¬ 
dial seed magnetic fields compatible with the currently avail¬ 
able upper limits leave detectable signatures on global prop¬ 
erties of the galaxy population, which can potentially be ex¬ 
ploited, in conjunction with the existing methods, to put 
tighter constraints on their intensities. 
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